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Yang-Mills theory is studied at finite temperature within the Hamiltonian approach in Coulomb 
gauge by means of the variational principle using a Gaussian type ansatz for the vacuum wave 
functional. Temperature is introduced by compactifying one spatial dimension. As a consequence 
the finite temperature behavior is encoded in the vacuum wave functional calculated on the spatial 
manifold x S^(I/) where L~^ is the temperature. The finite-temperature equations of motion are 
obtained by minimizing the vacuum energy density to two-loop order. We show analytically that 
these equations yield the correct zero-temperature limit while at infinite temperature they reduce 
to the equations of the 2-|-1-dimensional theory in accordance with dimensional reduction. The 
resulting propagators are compared to those obtained from the grand canonical ensemble where an 
additional ansatz for the density matrix is required. 

PACS numbers: 03.70.-|-k, 11.15.-q, 12.38.Aw, 

I. INTRODUCTION 

In recent years substantial efforts have been undertaken to study QCD non-perturbatively in the continuum. The 
main motivation is, of course, to understand the low-energy properties of QCD such as confinement and spontaneous 
breaking of chiral symmetry. Another motivation for these studies is the understanding of the phase diagram of QCD 
since the lattice approach can not deal with large chemical potentials due to the notorious sign problem. 

In Refs. [1-3] a variational approach to the Hamilton formulation of QCD in Coulomb gauge was developed.^ In 
the Yang-Mills sector this approach used initially Gaussian type variational ansatze for the vacuum wave functional. 
With the energy density calculated up to two-loops this variational approach has given a decent description of the 
infrared properties of the various propagators. Among other results one has found an infrared diverging gluon energy 
[1, 8] and a perimeter law for the spatial ’t Hooft loop [9] as well as a linear rising static color potential [8]. These 
are all signatures of confinement. Furthermore, it was shown that the inverse ghost form factor of the Hamiltonian 
approach in Coulomb gauge represents the color dielectric function of the Yang-Mills vacuum [10]. By the so-called 
horizon condition [11] this quantity vanishes in the infrared, which implies that the Yang-Mills vacuum behaves as a 
dual superconductor. 

The variational Hamiltonian approach in Coulomb gauge was extended in Refs. [12, 13] to finite temperatures 
by assuming a quasi-particle ansatz for the density operator of the grand canonical ensemble and minimizing the 
free energy. The obtained infrared exponents of the gluon and ghost propagators show a drastic change at the 
deconfinement phase transition from which a critical temperature of about Tq ~ 275... 290 MeV was obtained for 
SU(2) [13]. In Ref. [14, 15] the effective potential of the Polyakov loop was calculated within the Hamiltonian approach 
in Coulomb gauge. Thereby the temperature was introduced by compactifying one spatial dimension. Introducing the 
temperature into the Hamiltonian approach by compactifying a spatial dimension has the advantage that no ansatz for 
the density matrix (of the grand canonical ensemble) is required. The finite-temperature theory is fully encoded in the 
ground state on the spatial manifold IR^ x S^(L) with L~^ being the temperature, see Refs. [14, 15] for more details. 
Furthermore only this way of introducing the temperature makes the Polyakov loop accessible to the Hamiltonian 
approach, which uses the Weyl gauge Aq = 0. From the effective Potential of the Polyakov loop a critical temperature 
similar to that of the grand canonical ensemble is extracted. Furthermore, the correct order of the phase transition 
is obtained for both color groups SU(2) and SU(3). However, in the actual calculation of the Polyakov loop potential 
the propagators from the zero-temperature variational solutions were used, while, in principle, the finite-temperature 
self-consistent solutions should be used instead. Therefore in the present paper we extend the Hamiltonian approach 
to Yang-Mills theory in Coulomb gauge developed in Ref. [1] to the spatial manifold x S^(L). 

The rest of the paper is organized as follows: In the next section we derive the basic equations of a variational 
treatment of Yang-Mills theory in Coulomb gauge on to the spatial manifold x S^{L). Due to the explicit breaking 
of 0(3) symmetry on this spatial manifold this will require to generalize the ansatz for the vacuum wave functional 


^ The relation of the present approach [1] to previous variational calculations [4, 5] in Coulomb gauge is discussed in detail in Ref. [6]. 
Furthermore an analogous variational approach was recently developed in a covariant formulation in Landau gauge [7]. 
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of Ref. [1] distinguishing between the components of three vectors parallel and orthogonal to the compactified spatial 
dimension. In Sect. Ill we investigate the low- and high-temperature limits of this approach. The infrared analysis 
and the renormalization of our equations of motion are carried out in Sect. IV. In Sect. V we present our numerical 
results for the gluon and ghost propagators. Some concluding remarks are given in Sect. VI. 


II. FINITE-TEMPERATURE YANG-MILLS THEORY FROM A COMPACTIFIED SPATIAL 

DIMENSION 

Following Ref. [15] we introduce the temperature into the Hamiltonian approach to a quantum field theory by 
compactifying one spatial dimension, for which we choose the 3-axis. The partition function is then given 

Z(T) = lim Tre"'^«(^^ (1) 

l—¥00 

where Eq{L) is the energy of the vacuum state on the spatial manifold x S^(L) and I —>■ oo is the length of 
the uncompactified spatial dimensions. Note that in the present approach the whole temperature dependence is 
contained in the ground state energy Eq{L) and is completely encoded in the vacuum state. Excited states would 
only contribute when I is kept hnite, which corresponds to restricting one of the true spatial dimension to a finite 
length, as one encounters in the Casimir effect for two parallel infinite plates separated by a distance 1. 


A. Yang-Mills theory in Coulomb gauge on X S\L) 


On R^ X S^{L) the gauge field satisfies periodic boundary conditions in the third (compactified) dimension 

A(xi,T/2) = A(xj.,-T/2). (2) 

Here and in the following we denote by xj^ the component of x orthogonal to the compactified dimension. The 
Yang-Mills Hamiltonian in Coulomb gauge [16] 

dA = d±A± + (93^3 = 0 (3) 

on R^ X S^(L) is formally the same as on R^ 


= 2 d^x i^g^JX^UJAU + + He , 

except that the integration measure over the spatial manifold is defined by 

r r 


d a; := / d x_l 


dx3 . 


-L/2 


In Eq. (4) n^(a;) = —\5/5A‘l{x) is the momentum operator, B is the non-Abelian magnetic field and 


( 4 ) 

( 5 ) 


= Det(-Da) 


( 6 ) 


is the Faddeev-Popov determinant, where is the covariant derivative in the adjoint representation 

of the color group with being the structure constants. Finally, 

He = Y^d3a;^d^2/J^V“(x)J.4i^“''(x,y)p^(y) (7) 

is the Coulomb term, which arises from the elimination of the longitudinal components of the momentum operator 
by means of Gauss’ law. Here 




(x| (-D5)-i(-a^)(-D0)- 


|y) 


( 8 ) 


is the so-called Coulomb kernel, which is the non-Abelian generalization of the ordinary Coulomb law. 
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To find the vacuum state we use the variational principle with the trial vacuum wave functional [1] 

This ansatz preserves global color symmetry. The pre-exponential factor was chosen to cancel the Faddeev-Popov 
determinant in the functional integral for the expectation values 

{...)= j VAJAr[A]...^p[A]. (10) 

Due to the Coulomb gauge condition the variational kernel Wfc/(x, y) can be chosen transversal. On the spatial 
manifold IR^, due to 0(3) invariance, transversality implies 

Wfei(x,y) = tfe;(x)u;(x,y) , (11) 


where 

( 12 ) 

is the usual transversal projector in However, on x S^(L) the 0(3) symmetry is broken to 0(2) and the ansatz 
(11) has to be generalized to 

Wfci(x,y) = t^;(x)u;_L(x,y) -H4Ww||(x,y) . (13) 


Here 

4;(x) = (1 - 43 ) ( 4 * - (1 - Sis) (14) 

is the transversal projector in the subspace R^ orthogonal to the compactified 3-axis with Aj_ = df + 82 and 
furthermore 

4W=4;(x)-4W (15) 

is the corresponding orthogonal complement. There are now two independent variational kernels, W|| and a;_L, cor¬ 
responding to the compactified dimension and the orthogonal subspace, respectively. To simplify the notation we 
introduce an index Q G {||, T} to distinguish the compactified dimension from the space orthogonal to it. 

The projectors introduced above have the following properties 

4 w = 4m w. 4 w = 4m w > 

4fc(x) = 2, 4^ = 1. (16) 

With the trial ansatz (9), (13) the static gluon propagator is given by 

-Dfc?(x,y) = (4^(x)Hf(y)) = 4** (i:i^,(x,y) -f D^x.y)) (17) 


where 


Dki{^,y) = ^4(x)wQ^(x,y). 


- 


(18) 


By translational invariance the kernels a;Q(x, y) depend only on the relative coordinate x — y and it is convenient to 
work in momentum space, which is defined here by 


w(x) = / d^pe*P’"a;(p) := / d'^p±- ^ e' 


ip^Xj_ IWniCa 


Uj{pj_,UJn) , 


(19) 


n— — oo 


where d^p± = d'^p±/ (27r)^ is the integration measure in the momentum plane orthogonal to the 3-axis and = 27rn/L 
is the Matsubara frequency. From the representation (18) it is seen that wq(pj^,w„) with Q =|| and T, respectively 
are the energies of the transverse gluon modes parallel and orthogonal to the compactified spatial dimension. 
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FIG. 1. Ghost self-energy in rainbow-ladder approximation. The full and wavy line represent, respectively, the ghost and gluon 
propagator. 


B. Equations of motion 


The variational approach to Yang-Mills theory in Coulomb gauge on x S^(L) can be worked out in the same 
way as it was done in Ref. [1] on IR^ with the exception that the transverse degrees of freedom in are divided 
now into the degrees parallel and orthogonal to the compactified dimension, see Eq. (15). Accordingly there are two 
independent variational kernels wy and a;j_. In Ref. [13] the Coulomb term Hq (7) was shown to be irrelevant in the 
gluon sector and will therefore be ignored below. Following Ref. [1], but using the more general ansatz (13) for the 
vacuum wave functional, the variation of the energy, calculated up to two-loop order, results in the following two gap 
equations 

^|(p) = p^ + x|(p), OG{IU}, (20) 

where = p5_ + is the square of the 3-momentum 

p = Pj_ + a;„e3, = 27rn/T (21) 


and 

Xq(p) = y^cifci(p) / d^qqkqiG{p + q)G{q). (22) 

is the ghost loop, referred to in the present context as “curvature” [3]. Here G(p) is the ghost propagator defined by 

G = ((-Da)-i). (23) 


In momentum space the ghost propagator is conveniently expressed by its form factor d{p) defined by 


G(p) 


^(P) 

SP^ ' 


(24) 


The ghost form factor (i(p) expresses the deviation of QCD from QED and its inverse can be interpreted as the 
dielectric function of the Yang-Mills vacuum [10]. Calculating the ghost propagator with the vacuum wave functional 
(9) in the rainbow-ladder approximation results in the following Dyson-Schwinger equation (DSE) for the ghost form 
factor 


where the loop integrals 


d ^(P) = ^- (/j|(p)-b/^(p)) , 


I?{P) = Ne 



My(q)Pi d(p + q) 1 

p2 (p-bq)2 2 wq(p) 


(25) 


(26) 


arise from the ghost self-energy diagram, shown in fig. 1. DSEs, in general, are functional differential equations and 
their solutions are specified only after providing appropriate boundary conditions. For the ghost DSE we require (as 
on R^) the horizon condition, which on R^ x S^{L) reads 

d-\pj_ = 0,iOn = 0) = 0. (27) 


The two gap equations (20) and the ghost DSE (25) are coupled through the loop terms and have to be solved 
simultaneously. In the IR regime this can be done analytically following the by now standard IR analysis [17, 18], see 
Sect. III. The numerical solution for these equations in the whole momentum regime will be given in Sect. V. 
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III. ZERO- AND HIGH-TEMPERATURE LIMITS 

Before we present the numerical results let us investigate the above obtained equations of motion in the high and 
low temperature limit. 


A. The high-temperature limit 


In the high-temperature limit T = L~^ —>■ oo the Matsubara frequencies a;„ = 27rnT with n ^ 0 become infinite. 
Since for p —>■ oo all propagators vanish we need to keep in this limit in the above equations of motion only the lowest 
Matsubara frequency, ps = a;„ = 0, implying the replacement (cf. Eq. (19)) 




T 


d^ 


p± ■ 


For p 3 = 0 we have 


(pktli{q)pi) = 0 , 

^ y P3—q3—0 


(28) 


(29) 


so that the longitudinal part of the ghost self-energy vanishes, /j|(pj_,a;„ = 0) = 0, and the ghost DSE (25) reduces 
to 


d ^(P±,w„ = 0) 


= 0 ) 


(30) 


with 


^j'(P_L,W„ 


0) = N,T 


2 Pi*^A^±)pt d{p^ + (l^,Q) 1 

pi (Pu+Qu)^ 2u;i(q_L,0)' 


(31) 


This is precisely the ghost DSE in d = 2 spatial dimensions [19] provided we identify u;j_(qj_, 0) with the gluon energy 
and 


= d{pj_,0)VT 


(32) 


with the coupling constant and ghost form factor, respectively, in d = 2. As is well known the dimension of the 
Yang-Mills coupling constant g depends on the number of space-time dimensions. In d = 3 spatial dimensions g is 
dimensionless, while it has dimension y^energy in d = 2. Yang-Mills theory in d = 3 is scale free but the temperature 
does introduce a scale. Thus, a relation like Eq. (32) should have been expected. From Eq. (29) also follows that 
the “longitudinal” component of the curvature (22) vanishes in the high-temperature limit, X||(P±)‘^n = 0) = 0, and 
accordingly the longitudinal gap equation (20) reduces to 

Wy(Pj_,W„ = 0) = pi (33) 

implying that the gluon degree of freedom parallel to the compactified dimension becomes free at T —>- 00 . 

With Eq. (28) we find for the transversal curvature (22) 


X-l(p±,w„ = 0) 


g^TNc [ 2 d{pj_ + q^,0)d{q^,0) 

2 y qi d^(Pu + q^)" 


(34) 


which, with the identification (32), becomes the curvature in d = 2 spatial dimensions [19]. Accordingly the transversal 
gap equation (22) becomes that of two dimensions 

wi(Pi,Wri = 0) = pi-bxi(p,w„ = 0). (35) 


Thus for r —>■ 00 we find complete dimensional reduction: Yang-Mills theory in 3-1-1 dimension reduces in the high- 
temperature limit to the 2-|-l dimensional theory with an additional “Higgs field” given by the temporal component 
of the gauge field, Aq. In the present Hamiltonian approach this Higgs field does not show up since we work here in 
the Weyl gauge Aq = 0. 
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B. The zero-temperature limit 


As the temperature decreases more and more Matsubara frequencies have to be included. In the limit T —>■ 0 all 
Matsubara frequencies contribute and one has to carry out the whole infinite sum. For small T it is advantageous to 
carry out the summation over the Matsubara frequency by means of the relation [15] 


I OO 

2 


1 

27r 



OO 

dzf{z) 

k— — oc> 


(36) 


which follows from the Poisson resummation formula.^ 

In order to carry out the zero-temperature limit in the equations of motion given in section II consider a generic 
loop integral, which by means of Eq. (36) can be written as 


-j OO p p OO 

-p XI J = J d^pf{p±,P3) X (37) 


Here we have put z = ps and d^p = d^p/(27r)^ denotes the usual integration measure in 3-dimensional momentum 
space. As L —>■ OO the 0(3) invariance is restored and all propagators and correlation functions will depend only on 
the modulus of the 3-momentum, i.e. we have 

f{P±,P3) = f +pI^ = f{p) ■ (38) 

It is then convenient to use spherical coordinates in momentum space {p,6,(p). Carrying out the integration over 6 
we obtain from Eq. (37) 


1 

I 


X 

n— — oo 




pI 


■ OJt 


= / ^ 


k— — c 


sin(fcLp) 

kLp 


(39) 


Obviously the term fc = 0 reproduces the T = 0 result. The remaining terms vanish for L —>■ oo. With the relation 
(39) all equations given above reduce for A —>■ oo to the zero-temperature equations derived in Ref. [1]. 


IV. INFRARED ANALYSIS AND RENORMALIZATION 


In Ref. [18] an IR-analysis of the zero-temperature equations of motion was carried out using the power-law ansatze 
(for p 0) 

uj(p) ~ , d(p) ~ p~^ . (40) 

Assuming a bare ghost-gluon vertex and the horizon condition d~^(p = 0) = 0, i.e. /3 > 0, one obtains from the ghost 
DSE the sum rule 


a = 2P + 2-d, (41) 

where d is the number of spatial dimensions. Erom the gap equation one finds then two solutions for d = 3 

13=1, 0.7952 (42) 

and a single solution for d = 2 

/3 = 0.4 . (43) 

The IR exponents extracted from the numerical (T = 0) solutions obtained in Refs. [1, 8] (see also Ref. [13] and [19]) 
are in excellent agreement with the analytic results of this infrared analysis. Erom the result of the previous section it 


^ When Eq. (36) is applied to functions f{uin) known only numerically for the discrete Matsubara frequencies u)n it is understood that 
the function f{z) arises from f{io„) by smoothly interpolating the latter in the intervals between neighboring Matsubara frequencies. 
In fact, for L —> oo, where the Poisson resummation (36) is convenient, the spacing between neighboring Matsubara frequencies tends 
to zero. 
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is clear that the propagators on x S^(L) with the Matsubara frequency n = 0 will exhibit the d = 3 IR exponents 
(42) in the low temperature regime and the d = 2 IR exponents (43) in the high-temperature regime. In this context 
it is worth mentioning that in the finite-temperature extension of the Hamiltonian approach to Yang-Mills theory 
in Coulomb gauge within the usual canonical ensemble given in Ref. [13] the IR-analysis can be carried out in the 
high-temperature limit yielding an IR exponent of /3 = 0.5. In the corresponding numerical studies the deconfinement 
phase transition was accompanied by a rapid change of the IR exponents from the (T = 0, d = 3)-values (42) to the 
high-temperature value of P = 0.5, close to the 2-dimensional value (43). 

As shown in the previous section for L —>■ oo (T —>■ 0) the equations of motion of the present approach reduce 
indeed to the zero-temperature equations of Ref. [1]. By compactifying one spatial dimension no (additional) UV 
singularity can be introduced. Furthermore, for large momenta p ^ T the temperature becomes irrelevant. Thus 
our equations of motion on x S^(L) will have the same UV behavior as on R^ and can be renormalized as in the 
zero-temperature case. For T = 0 the ghost loop x{p) (22) is linearly UV divergent, while the loop integrals I^{p) 
(26) in the ghost DSE are logarithmically divergent. To renormalize the ghost DSE (25) we use Eq. (15) to rewrite 
it in the form 


d ^(p) = - - Vc [ d^q [l - (pq)^] 
9 J L 


djp + q) 1 

(P + q)^ 2u;j^(q) 


-Nr 


3 Mj-(q)Pj d(p-bq) 


d-^q 


(p + q)^ V2w|i(q) 


■ (44) 


To renormalize the ghost DSE it would, in principle, be sufficient to subtract from the first loop integral Eq. (44) its 
zero-temperature counterpart at a renormalization scale pd- Although being correct this procedure is not feasible in 
an iterative solution of the ghost DSE. Therefore we renormalize the ghost DSE by subtracting the first (temperature- 
dependent) loop integral at a scale pd, which is chosen in the UV such that the horizon condition (27) is fulfilled.^ 
Analogously we renormalize the curvature x'^(p) by subtracting it at a scale p^ 

X‘5(p)^x'^(p)-X'^(dx)- ( 45 ) 

We choose the vectors Pd = (Pd,_L) dd,™) and = {Px,A.i 9x,n) for given subtraction points pd = \Pd\ and p^ = |p^| 
to be parallel to the external momentum p, i.e. we put 

Pd = ddP, Px = dxP> 

For a more general discussion of the renormalization of the Hamiltonian approach in Coulomb gauge see Ref. [20], 
and Ref. [13] for the renormalization at finite temperature. 


V. NUMERICAL RESULTS 

Our equations of motion (20) and (44) represent a set of coupled integral equations with respect to the transverse 
momentum pj^. Each of these integral equations labeled by (the index of) the Matsubara frequencies a;„ is solved 
using standard methods, see e.q. Ref. [20]. The equations belonging to different Matsubara frequencies are coupled 
through the loop terms. Eor small temperatures many Matsubara frequencies need to be included and the number of 
integral equations becomes large. To reduce the numerical expense we solve the integral equations explicitly only up 
to a maximum Matsubara frequency , while for > oJumax tbe extrapolation 

d(p±, w„) = (47) 

is used. This relation is exact for 0(3) invariant solutions, which depend on pj_ and only through the combination 
\/p'j_ + uJn- A-S expected and as a test calculation with n = 20 Matsubara frequencies confirms, our solutions become 
quasi 0(3) invariant for momenta large compared to the temperature, i.e. in particular for 3> T or 27rn 1. In 
fact this test calculation shows that the interpolation (47) works very well for Umax > 4. In our numerical calculations 
we have chosen rimax = 4 and included in the sums up to n = 100 Matsubara frequencies. 


® This introduces de facto a temperature-dependent renormalization constant. The same renormalization procedure was used in the study 
of the grand canonical ensemble in Ref. [13]. 





(a) (b) 

FIG. 2. The ghost form factor d{p±,u}n) (a) as function of the transverse momentum p± for the lowest Matsubara frequency 
n — 0 and for various temperatures decreasing from bottom to top, and (b) at fixed temperature for the first Matsubara 
frequencies n — 0,... ,3. 



FIG. 3. (a) the transversal Ljj±{p±,oJn.) and (b) the longitudinal cu||(p^,a;„) gluon energies for the lowest Matsubara frequency 
n = 0 as function of the transversal momentum for various temperatures decreasing from bottom to top. 


In previous studies [1, 8] at zero temperature fj-d = 0 was chosen which allows to explicitly implement the horizon 
conditon d“^(0) = 0 into the renormalized ghost DSE. Depending on the starting function of the iteration one obtains 
then either the /3 = I or /3 si 0.8 solution. Here we follow a different road, which actually turns out to be more 
efficient: We choose pd in the ultraviolet and fine tune d~^{p,d) such that the horizon condition d“^(0) = 0 is fulfilled. 
For given ^ 0 there are two distinct values for d~^{p,d) for which the horizon condition is realized, one resulting 
in the (3 = 1 solution, the other in the /3 ~ 0.8 solution. We consider here only the /3 = 1 solution, which at T = 0 
yields a linearly rising non-Abelian Coulomb potential. 

Fig. 2a shows the ghost form factor for the Matsubara frequency n = 0 as function of the “transverse” momentum 
pj_ = Ipj^I and for various temperatures. In the high-temperature limit d{p±,uJn = 0) in fact becomes the T = 0 
solution of d = 2 spatial dimensions and acquires in particular the IR exponent (43) of that solution. As the 
temperature is lowered in the mid-to-small-momentum regime a second power-law emerges with an exponent given by 
the IR exponent (42) /? = I of the T = 0 solution of d = 3. This power-law extends to smaller and smaller momenta 
as the temperature is lowered, while the power law of the (T = 0, d = 2) solution stills persists in the deep IR for any 
finite temperature and disappears only in the zero-temperature limit (L —> oo). 
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FIG. 4. (a) the transversal ui±{jp±,ujn) and (b) the longitudinal uj\\{p±,uin) gluon energies at finite temperature as function of 
the transversal momentum and for the first Matsubara frequencies n = 0,..., 3. 


An analogous temperature behavior is found for the transversal and longitudinal gluon energies shown in Fig. 3. 
At high-temperatures u;_L(p_L,a;„ = 0) approaches the {T — 0, d = 2)-solution, while W||(pxjW„ = 0) is IR finite. As 
the temperature is lowered, in the mid-to-small momentum regime for both ui±(jj±,ujn = 0) and uj\^(j)±,ujn = 0) a 
power-law behavior uiQ{p±,uJn = 0) ^ with the d = 3 (T = 0) IR exponent a = 1 emerges. However, in the deep 
IR ujQ{p±,uin = 0) behaves still as in the high-temperature limit, i.e. u}±{p±,uJn = 0) shows the power law with the 
d = 2 IR exponent a = 0.8 while uj^\{p±,ujn = 0) goes to a plateau for p± —>• 0. With decreasing temperatures the 
d = 3-power-law covers an increasing portion of the IR-momentum regime and for T —>■ 0 both uj±{p±,u!n = 0) and 
W|| {p±,ujn = 0) merge to the 0(3) invariant (IR-diverging) (T = 0, d = 3)-solution, in accord with our findings for the 
zero-temperature limit given in Sect. III. Finally Fig. 4 shows the gluon energies uj± and W|| at finite temperatures 
and for the first few Matsubara frequencies n = 0,..., 3. For n > 0 both gluon energies as well as the ghost form 
factor [Fig. (2b)] are infrared finite. 


VI. CONCLUSIONS 


In this paper we have studied finite-temperature Yang-Mills theory by compactifying a spatial dimension to a 
circle S^(L) and interpreting its circumference L as the inverse temperature. In this approach the finite-temperature 
field theory is entirely encoded in its vacuum state calculated on spatial manifold x S^(L). The vacuum wave 
functional on the x S^(L) was determined by the variational principle using a Gaussian type trial state which 
was generalized to have different variational kernels for the gauge field degrees of freedom parallel and orthogonal 
to the compactified dimension. The horizon condition was imposed for the zeroth Matsubara frequency. At high 
temperatures T —)> oo the self-consistent solutions approach those of the d = 2 dimensional theory at T = 0, as 
expected from dimensional reduction. As the temperature is lowered the self-consistent solution exhibits in the 
medium-to-small momentum regime a power-law behavior with the IR exponent of the (T = 0,d = 3)-solutions, while 
in the deep IR the (T = 0,d = 2)-power-law still persists. The IR momentum range of the latter decreases with 
decreasing temperature and vanishes only in the zero-temperature limit. 

The present approach to finite-temperature quantum field theory is advantageous compared to the standard ap¬ 
proach based on the grand canonical ensemble since it requires no ansatz for the density matrix. The presently 
obtained self-consistent variational solution of Yang-Mills theory on R^ x S^{L) are required as input for a fully self- 
consistent calculation of the Polyakov loop at finite temperature in the Hamiltonian approach. This will be subject 
to future work. 
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